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ABSTRACT 

Our PACS and SPIRE images of the Aquila Rift and part of the Polaris Flare regions, taken during the science demonstration phase of Herschel 
discovered fascinating, omnipresent filamentary structures that appear to be physically related to compact cores. We briefly describe a new multi- 
scale, multi-wavelength source extraction method used to detect objects and measure their parameters in our Herschel images. All of the extracted 
starless cores (541 in Aquila and 302 in Polaris) appear to form in the long and very narrow filaments. With its combination of the far-IR resolution 
and sensitivity, Herschel directly reveals the filaments in which the dense cores are embedded; the filaments are resolved and have deconvolved 
widths of ~ 35" in Aquila and ~ 59" in Polaris (~ 9000 AU in both regions). Our first results of observations with Herschel enable us to suggest 
that in general dense cores may originate in a process of fragmentation of complex networks of long, thin filaments, likely formed as a result of an 
interplay between gravity, interstellar turbulence, and magnetic fields. To unravel the roles of the processes, one has to obtain additional kinematic 
and polarization information; these follow-up observations are planned. 
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1. Introduction 



Herschel (Pilbratt et al. 2010) provides the best opportunity to 
study the earliest stages of star formation and the origin of the 
stellar initial mass function (IMF). Prestellar cores and young 
(Class 0) protostars emit the bulk of their luminosities at wave- 
lengths 80-400 /i m, which makes th e Herschel imag ing in- 
strum ents PACS (IPoglitschetal .11201 Oh and SPIRE (Griffi n et all 
2010) perfect for performing a census of these objects down to 
0.01-0. IMq in the nearby (distances D ;£ 500 pc) molecular 
clou d complexes. Our Herschel key project (Gould Belt sur- 
vey, lAndre et al]|2010h aims at probing the link between diffuse 
cirrus-like structures and compact cores with the main scientific 
goal to elucidate the physical mechanisms of the formation of 
prestellar cores out of the diffuse medium, which is crucial for 
understanding the origin of stellar masses. 

During the science demonstration phase of Herschel obser- 
vations we imaged two large fields, the Aquila Rift and part of 
the Polaris Flare (hereafter called Aquila and Polaris fields for 
brevity), with the aim of better understanding star formation in 
extremely different environments. The high-quality far-infrared 
images are extremely rich in new information on the earliest 



phases of star formation, which will be explored in full detail 
in the future. In this short paper we focus on merely one dis- 
covery made by Herschel: the impressive, ubiquitous network 
of filamentary structures in the interstellar medium (ISM) and 
its remarkably close physical relationship with the objects that 
appear to form within the filaments. 

2. Observations and data reduction 
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The PACS instrument is described by Poglit sch et alj d2010l) ; the 
SPIRE instrument, its in-orbi t performance, and its scientific ca- 
pabilities are described by Griffin et al] d2010h . and the SPIRE 
astronomical c alibration methods and accuracy are outlined by 
ISwinvard et alJ (l2010t) . Our Herschel o bserv ations of the Aquila 
field are presented bv iKonvves et al.1 (l2010h : we describe here 
our observations of the Polaris field. 

The entire Polaris field (8 square degrees) was observed 
on 2009 October 23, at 70, 160, 250, 350, and 500 /an in the 
PACS/SPIRE parallel mode at 60 "s -1 , performing cross-linked 
scans in two orthogonal directions. The PACS raw data were re- 
duced with the HIPE 3.0.455 software with standard steps of the 
default pipeline. The baseline was subtracted by high-pass filter- 
ing with a ~1° median filter width (the full leg length was 2?5). 
Multi-resolution median transform deglitching and second order 
deglitching were applied, and final maps were computed by the 
HIPE's MADmap method. The SPIRE parallel-mode data were 
reduced with HIPE 2.0 using pipeline scripts delivered with that 
version, modified to include data that were taken during the tele- 
scope's turnaround at the scan ends. A median baseline was ap- 
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plied to the maps on individual scan legs and a "naive mapper" 
was used to create images. The images were projected to the 
nominal 3'.'2, 3'.'2, and 6'.'4 pixels for PACS and 6", 10", and 14" 
pixels for SPIRE wavelengths. In our analysis, we adopted cir- 
cular Gaussian point-spread functions (PSFs) of 8'.'4, 9'.'4, 13'.'5, 
1871, 2572, and 3679 at the wavelengths (although the beams 
are actually somewhat elongated). Before extracting objects, all 
images were resam pled to the same 6" pixel size using SWarp 
(Berli n et al.ll2002l) . 



3. Source extraction method 

Here, we summarize the multi-scale, multi-wavelength source 
extraction algorithm that has been developed and extensively 
tested by A. M. at CEA Saclay. The main motivation for the de- 
velopment was the absence of any multi-wavelength extraction 
technique for use in our Herschel projects. Full description of 
the method and code (called getsources) as well as of the bench- 
marking results for this and several other extraction technique^] 
will be given elsewhere (Men'shc hikov et aDl2010h . 

The main idea of the method used in this work is to analyze 
decompositions of original images (at each wavelength) across 
a wide range of spatial scales separated by only a small amount 
(typically ~ 5%). Each of those "single scales" are cleaned of 
noise and background by iterating to appropriate cut-off levels, 
then re-normalized and summed up over all wavelengths in a 
combined single-scale detection image. The main advantage of 
this algorithm is in its multi-wavelength design: the same com- 
bined detection image across all wavelengths eliminates the need 
of matching multiple catalogs obtained with different angular 
resolutions and reduces associated measurement errors. Besides, 
fine spatial decomposition improves the detection of even the 
faintest objects and aids in better deblending of crowded regions. 

The decomposition is done by convolving the original im- 
ages with circular Gaussians and subtracting them from one an- 
other (this can be called successive unsharp masking): 



Ij(A) = Gj-i * 1(A) - Gj * 1(A) (j = 1,2, ...,A0, 



(1) 



where 1(A) is the original image at a wavelength A, Ij(A) are 
"single-scale" decompositions, Gj are the smoothing Gaussians 
(Go is a two-dimensional delta function) with PSF $ FWHM $ 
image size, FWH1VL; = f s x FWHM ; -_i, / s > 1 is a scale factor 
(usually f s ~l .05), PSF is the smallest beam over all A, and the 
number of scales N depends on the value of f s . 

To separate the signals of objects from the noise and back- 
ground contributions, we iterate to 5<Tj(A) cut-off levels in each 
of Ij(A) (the standard deviation ctj(A) is computed outside the 
objects) and zero out all pixels with intensities below that level, 
producing clean images I C XA). An advantage of this single-scale 
cleaning is that the noise contribution depends very significantly 
on the scale (e.g., at large scales the small-scale noise is di- 
luted, whereas large objects are best visible). In effect, in a full 
clean image I C (A) =£j I C :(A) one would see structures deeper 
than in 1(A). The above algorithm is applied to images at each 
band separately and the clean images are used to produce com- 
bined single-scale detection images l z -Yj\ f X I%A) after nor- 



malization to a similar maximum intensity^ Objects in the com- 
bined detection image are identified by tracking their appearance 
and "evolution" from small to large scales, and their footprints 
are determined. The objects are background-subtracted and de- 
blended, and their sizes and fluxes are measured in the original 
observed images at each wavelength. For more details, we refer 
interested readers to Appendix A. 

4. Results and discussion 

To derive properties of the compact objects in the Aquila and 
Polaris fields, we performed source extraction with getsources 
(Sect. [3]) and made a careful analysis of the resulting extraction 
catalogs, selecting most relia ble objects and determining core 
mass function for both fields dKonvves et alJl2010l : lAndre et all 
2010). To assess completeness levels of our observations, we 
created the synthetic skies with the actual observed background 
of our images (after removing all extracted objects) and pop- 
ulated them wit h cores and protostars from spherical radiative 
transfer models dMen'shchikov et al]|2010h . then performed ex- 
tractions on the skies. 

The filamentary structures are clearly visible in the left pan- 
els of Fig. Q] where both the Aquila and Polaris fields are shown 
as 3-color composites created from the images at 500, 160, and 
70 fj.m. With red color, the images highlight strong density en- 
hancements and low temperatures, whereas the material within 
two bright hot regions (HII region W 40 in the north and around 
MWC 297 and Sh 2-62 stellar group in the south) of the Aquila 
field is visible in green, blue, or white. From this one could im- 
mediately conclude where to look for prestellar cores, and it is 
indeed the red filaments where the reddest cores are visible. The 
colors are not very pronounced in the Polaris field, suggesting 
much smaller temperature and density variations across the field, 
even in its densest filaments. To better visualize and characterize 
the filaments, we utilized several approaches. 

High-contrast images of the cores and filaments in both fields 
were obtained as part of our source extraction technique by ap- 
plying the successive unsharp masking (Eq. Q} with / s =1.05, 
Gi=9739, and N-99. In particular, the right panels in Fig. Q] 
show that the visibility of the very narrow and long filaments 
is dramatically improved in some of our "single-scale" images 
(which effectively select a narrow range of scales, with both 
larger and smaller scales filtered out). For the sake of illustration, 
a 40" scale (in which the structures are best visible) is displayed 
in Fig. [1] The compact objects are populating the filaments al- 
most like pearls on threads in a necklace, giving them a nodu- 
lar, lumpy appearance. The same decomposition at all Herschel 
wavelengths proved that the same structures are seen across the 
entire range, from 70 to 500 yum. There are, however, noticeable 
changes in the appearance of the filaments in Aquila (around 
W40) at the shortest wavelengths (Fig.|2j upper left), presum- 
ably due to increasing opacity and to changes in temperature of 
some of the filaments (complex three-dimensional structures). 

We also extracted filamentary structure from the observed 
imag es with the M CA software (morphological component anal- 
ysis, IStarck et al]|200 4). The idea of the method is to morpho- 
logically decompose a signal into its building blocks, which 
can be represented by isotropic wavelets, ridgelets, or curvelets. 



1 The Gould Belt consortium is testing several source extraction al- 
gorithms, such as ga ussclumps (Stutzki & Guesten 1990), clumpfind 
i Williams et al. 1994), sextractor (Bertin & Arnouts 1996), derivatives 
i Molinari et al. 2010), and a few others (see Appendix A for a very brief 
summary), using simulated skies of various degrees of complexity. 



2 At this detection step, it is neither possible nor necessary to prop- 
erly preserve the spectral behavior of sources at single scales in the 
combined detection image, as the processing is done at the entire im- 
age level before finding objects; spectral properties of the latter will be 
determined from the original images at the measurement step. 
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Filamentary structures were separated from the more isotropic 
features by applying both the curvelet and wavelet transforms 
in 100 iterations (with exponentially decreasing thresholds) to 
achieve convergence in both the wavelet and curvelet represen- 
tations of the original images . The curv elet component (Fig. [2] 
right panels; see also I Andre et al.ll2010h shows filaments very 
similar to those displayed in Fig. Q] and left panels of Fig. |2j we 
also get analogous results wh en processing the col umn density 
images (the latter are shown in lKonvves et al.ll2.01Qt) . 

Measurements of the thickness of several high-contrast fila- 
ments in our column density maps (obtained with a uniform 39" 
resolution of the 500 //m band) give 48" for Aquila and 66" for 
Polaris (FWHM), uncertain to within 12"; when deconvolved, 
they are 28" and 53". If we adopt the distances to th e regions of 
260 pc dStraizvs et alJl2003b iBontemps et alJl20lol) and 150pc 
(Heithausen & Thaddeus 1990), the deconvolved linear widths 
are ~ 7300 AU and 8000 AU. Measurements of the widths at 
several locations across a few well-behaved filaments (those not 
blended with the other filaments) in the observed images re- 
sulted in the values of 30, 36, 42, 42, 54" (for Aquila at the 
PACS+SPIRE wavelengths) and 63, 64, 69" (for Polaris at the 
SPIRE wavelengths). The filaments are more resolved at PACS 
wavelengths and the values would imply average deconvolved 
widths of 35" for Aquila and 59" for Polaris (9000 AU in both 
fields). The largest FWHM sizes of the extracted cores are com- 
parable to the filaments' widths (within ~50%) for both fields, 
confirming the visual impression that the objects are largely con- 
fined to the filaments they are embedded in. There seems to be 
a wide distribution of the lengths of the filaments; some of them 
are as long as ~(f.5 (Figs. Q~l |2|, which would correspond to a 
few pc. It is clear, however, that one needs to perform a more 
systematic and complete study of the filaments' properties; that 
will be a subject of a future work. 

The Aquila and Polaris fields present very different physical 
conditions; the former shows very active star formation, while 
the latter displays almost none. Indeed, we find many extremely 
cold and dense filaments (T d « 7.5-15 K, N Hl ~ 5 x 10 20 - 
1.4 x 10 23 crrT 2 ) in Aquila; most of them are gravitationally un- 
stable and fragme nted into several hundreds of prestellar cores 
(Andre et al. 2010). In contrast, the Polaris field contains warmer 
filaments with much lower column densities (Ta « 10-15 K, 
N H2 ~ 3 x 10 20 -8.6 x 10 21 cm -2 ) and w e find no clear examples 
of pr estellar cores (Andre et al. 20 1(J IWard-Thompson et alj 
2010). To illustrate the very close relationship between the fil- 
aments and detected objects, we overplot locations of the latter 
on the curvelet transform for both fields. An inspection of Fig. [2] 
shows that all objects detected in Aquila and Polaris are indeed 
located within very narrow filaments. Our maps of a stability pa- 
rameter of the filaments (their mass per unit length) suggest that 
the Polaris filaments are stable, w hile many dense fil aments in 
Aquila are gravitationally unstable dAndre et al.l F2010). 

Several possible models for the formation of filamentary 
cloud structures have been proposed in the literature. In par- 
ticular, numerical simulations of supersonic magnetohydrody- 
namic turbulence in weakly magnetized clouds always generate 
complex systems of shocks, which fragment the gas into high- 
density sheets, filaments, and cores (e.g. Pad oan et alj 12001). 
Filamentary structures are also produced in turbulent simula- 
tions of more strongly magnetized m olecular clouds, includ - 
ing the effects of ambipolar diffusion (Nakam ura & Lil 12008). 
Another possibility is that the filaments originate from grav- 
itational fragmentation of sheet-like cloud layers (Nagaietal. 
1998). Interestingly, filaments can result from purely gravita- 
tional amplification of a random field of initial density fluc- 



tuations as is the case in large cosmological simulations (e.g. 
lOcvirk et alJl20 08). Our measurements of the radial profiles of 
selected filaments in Aquila and Polaris suggest that magnetic 
fields may play an important role in shaping the structure of the 
observed filaments, if not in forming the filaments. Indeed, they 
indicate that the density falls off approximately as p oc r~ L5 (in 
Aquila) or p oc r~ 2 (in Polaris) away from the axes of the fil- 
aments. In both cases, the measured radial density profiles are 
much shallower than the steep p oc r ~ 4 profile of unmagnetized 
hydrostatic filame nts fostrikerlll96 4). as already pointed out by 
lLadaetalJ d 1999b and lJohnsto ne & Bally ( 1999) for two specific 
filaments. Models of unmagnetized equilibrium filaments are 
therefore inconsistent with observations. By contrast, models of 
equilibrium filaments with primarily toroidal or helical magnetic 
fields can account for radial profiles ra nging from p oc r~ l t o 



p oc r~ 2 , in agreement with observations (Fiege & Pudritz 200' 

Filaments are seen in numerous astronomical images and the 
filamentary stru cture of molecular clouds has been known for 
some time (e.g. [Schneider & Elmegreenl l 19791; iGoldsmith et al.1 
2008). However, Herschel with its combination of the far-IR res- 
olution and sensitivity directly reveals the filaments in which the 
cores are forming. Our first observations demonstrate the fas- 
cinating, omnipresent filamentary density structure of the cold 
ISM and its extraordinarily intimate physical relationship with 
the objects that form within the filaments. They enable us to 
suggest that in general dense cores may originate in a process 
of fragmentation of long, thin filaments, presumably formed as 
a result of an interplay between gravity, int erstellar turbulence , 
and magnetic fields (see also discussion in lAndre et al.l 12010). 
To unravel the role and relative importance of these processes, 
one needs additiona l kinematic and polarization information (cf. 
Matt hews et al.l200lh : these follow-up observations are planned. 
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Fig. 1. Composite 3-color images (fe/f) produced from the observed images of the Aquila (fop) and Polaris (bottom) fields (panels 
are 3!3x3?3, or ~15 pc and ~9 pc at their distances of 260 pc and 150 pc, respectively; see Sect. 3); the red color (SPIRE band 
at 350/im), highlights colder areas and objects, whereas the green and blue colors (PACS bands at 160 /im and 70/im) show 
regions with progressively hotter radiation fields. High-contrast "single-scale" decompositions (right) of sub-fields in Aquila (top) 
and Polaris (bottom) combined from those in all SPIRE bands for better visibility, display intensity distribution (MJy/sr) within a 
narrow range of spatial scales around 40"; the sub-fields of 1°2 x 1°2 are outlined by white squares in the left panels. Clearly visible 
is a very tight association of the narrow filaments with compact cores. 
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Appendix A: Extraction techniques 

A.1. Existing source extraction algorithms 

Here we summarize (very briefly) the concepts of different tech- 
niques, to place getsources described in Sect.[3]in a wider con- 
text. The algorithms trying to solve the same problem of source 
extraction originated from different ideas. Note that they have 
also been developed (oriented) for use in different areas of as- 
tronomy, thus their performance for a specific project must be 
carefully tested before an appropriate method can be chosen. 

Stutzki & Guesten (1990)'s gaussclumps (originally created 
for position-velocity cubes) fits a Gaussian profile to the bright- 
est peak, subtracting the fit from the image, then fitting a new 
profile to the brightest peak in the ima ge of residuals, iteratin g 
until some termination criteria are met. William s et al.l (fl994)'s 
clumpfind contours an image at a number of levels, starting 
from the brightest peak in the image and descending down to 
a minimum contour level, marking as clumps along the way 
all connected areas of p ixels that are above the contour level. 
Berti n & Arnoutsl dl996l) 's sextractor estimates and subtracts 
background, then uses thresholding to find objects, deblends 
them if they overlap, and measures their positions and sizes 
using intensity moments. CUPIEqs reinhold identifies pixels 
within the image which mark the edges of clumps of emission, 
producing a set of rings around the clumps. After cleaning noise 
effects on the edges, all pixels within each ring are assumed to 
belong to a single clump. CUPID's fellwalker ascends image 
peaks by following the line of the steepest ascent, considering 
every pixel in the image as a starting point for a walk to a signif- 
icant peak , marking along the way all visited pixels with a clump 
identifier. Motte et al. (2007)'s mre-gcl combines cloud filter- 
ing t echniques based on wa velet multi-resolut ion algorithms 



I e.g., Starck & M urtagh 2006) with saussclumps. M olinari et al] 
2010)'s derivatives analyzes multi-directional second deriva- 
tives of the original image and performs curvature threshold- 
ing to isolate compact objects, then fits variable-size elliptical 
Gaussians (adding also a planar background) at their positions. 
Another method that defines cores in terms of connected pixels 
is csar, which was developed for use with BLAST and Herschel 
(Harry et al. 2010, in preparation). 

Whereas clumpfind, reinhold, fellwalker, and csar merely 
partition the image between objects not allowing them to over- 
lap, gaussclumps, sextractor, and mre-gcl can deblend overlap- 
ping objects, which is quite essential for obtaining correct results 
in crowded regions. None of the methods was designed to handle 
multi-wavelength data, making it necessary to match the cata- 
logs obtained at different wavelegths using an association radius 
as a parameter. 



to the object across all scales. The segmentation masks expand 
toward larger scales, and the evolution of each object's mask is 
followed, as is appearance of new objects at any scale and dis- 
appearance of those which become too faint at the current and 
larger scales. When two or more objects touch each other in a 
single-scale image, the segmentation masks are not allowed to 
overlap, but overlapping does happen between objects of differ- 
ent scales. The largest extent of any source defines its footprint, 
and this is determined at the scale where the object's contrast 
above the cut-off level is at maximum. The scale itself provides 
an initial estimate for the object's FWHM size. 

The positions of sources are computed from the first mo- 
ments of intensities in a combined detection image at a range of 
single scales, from where an object first appeared and to the scale 
twice as large. The objects' sizes are computed from the first and 
second intensity moments in the original background- subtracted 
image. The background subtraction is done by linearly interpo- 
lating pixel intensities off the observed image under the foot- 
prints, in the four main directions (two axes and two diagonals), 
based on the pixel values just outside the footprints. Our itera- 
tive deblending algorithm employs two-dimensional shapes with 
peak intensities and sizes of the extracted objects in order to di- 
vide the intensity of a pixel between surrounding objects accord- 
ing to the fraction of the shapes' intensities at the pixel. For the 
shapes we adopted a two-dimensional analog of the Gaussian- 
like function / (1 + (r/r Q ) 2 )- a (Moffat 1969) with a = 10. 

The end result of the processing is an extraction catalog (one 
line per object) containing coordinates of all detections (inde- 
pendent of A) and estimates of the objects' S/N ratios, peak and 
total fluxes (with their uncertainties), and sizes and orientations 
for each wavelength. In addition, getsources produces catalogs 
of all possible colors, as well as the azimuthally-averaged inten- 
sity profiles (their full, background-subtracted, and deblended 
versions) and deblended images for each object. 



A.2. More details on the new method 

In getsources the extraction of objects is performed in each of 
the combined detection images by going from the smallest to 
the largest scales and finding segmentation masks of the objects 
at each scale using the tint fill algorithm (Smith 19790- The 
masks are the areas of connected pixels in a segmentation im- 
age, and the algorithm fills the pixels' values with the number 
of a detected object and allows tracking of all pixels belonging 



3 CUPID is a source extraction software package developed by the 
STARLINK team for use with the SCUBA2 surveys; it is a general 
wrapper to which additional methods can be added. See documentation: 
http://docs.jach.hawaii.edu/star/sun255.htx/sun255.html 

4 Available at http://portal.acm.org/citation.cfm?id=800249.807456 
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Fig. 2. Composite 3-color images left of the sub-fields of Aquila (top) and Polaris (bottom) produced from the high-contrast "single- 
scale" decompositions (red comes from all SPIRE bands, green and blue correspond to the PACS bands at 160 and 70 pm, re- 
spectively); the sub-fields and decompositions are the same as in Fig. Q] Curvelet components (right) extracted from the original 
SPIRE 350/im images overlaid with ellipses for selected cores extracted by getsources (shown are starless cores: those detected 
at S/N>7.5 levels in at least two SPIRE bands and not detected in emission or detected in absorption in the PACS 70pm band, 
with SED dust temperatures T& < 18 K; see also lKonvves et al . 2010; Andre et al. 2010). Most objects lie in the narrow and long 
filaments spanning orders of magnitude in intensity (MJy/sr). 



